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Abstract Methanogenic archaea use a [NiFe]-hydrogenase, Frh, for oxidation/reduction of F420, 
an important hydride carrier in the methanogenesis pathway from H2 and CO2. Frh accounts for 
about 1% of the cytoplasmic protein and forms a huge complex consisting of FrhABG heterotrimers 
with each a [NiFe] center, four Fe-S clusters and an FAD. Here, we report the structure determined 
by near-atomic resolution cryo-EM of Frh with and without bound substrate F420. The polypeptide 
chains of FrhB, for which there was no homolog, was traced de novo from the EM map. The 
1 .2-MDa complex contains 1 2 copies of the heterotrimer, which unexpectedly form a spherical 
protein shell with a hollow core. The cryo-EM map reveals strong electron density of the chains of 
metal clusters running parallel to the protein shell, and the F42o-binding site is located at the end of 
the chain near the outside of the spherical structure. 
DOI: 10.7554/eLife.0021 8.001 



Introduction 

[NiFe]-hydrogenases are microbial enzymes that heterolytically cleave H2, after which the electrons 
from a hydride ion are reversibly transferred to electron carriers. These enzymes are involved in many 
metabolic pathways in microbial ecosystems, notably in methanogenesis. The members of the [NiFe]- 
hydrogenase family are composed of at least two subunits, large (~60 kDa) and small (~30 kDa), which 
contain a [NiFe] dinuclear metal center and three iron-sulfur (Fe-S) clusters, respectively. The family is 
classified into five distantly related groups based on a phylogenetic analysis of their primary structures 
{Vignais and Billoud, 2007; Constant et al., 2011). The cytoplasmic F42o-reducing [NiFe]-hydrogenase 
(Frh) belongs to the group 3 [NiFe]-hydrogenases, which catalyze the reversible reduction of the soluble 
hydride carrier, NAD(P) or F420. F420 is a deazaflavin derivative that acts as an important hydride accep- 
tor/donor in the central methanogenic pathway. The group 3 hydrogenases contain an additional 
subunit that interacts with the soluble coenzymes, for example, the iron-sulfur flavoprotein FrhB. To 
date, the crystal structures of [NiFe]-hydrogenases only from group 1 have been determined: six derived 
from sulfate-reducing bacteria, including two orthologs harboring a [NiFeSe]-center as active site 
(Vb/faeda et a/., 1995; Higuchi et al., 1997; Montet et al.. 1997; Garcin et al., 1999; Matias et al.. 2001; 
Marques et al.. 2010), one from a photosynthetic bacterium {Ogata et al.. 2010), and three from 
hydrogen-oxidizing bacteria and Escherichia coli, which in contrast to the other species are oxygen- 
tolerant and have an usual [4Fe-3S] proximal cluster {Fritsch et al., 2011; Shomura et al.. 2011; 
Volbeda et al., 2012). 

Frh is a key enzyme in the metabolism of methanogenic archaea {Thauer et al.. 2010). The reduction 
of carbon dioxide to methane is an overall eight-electron reduction process involving the oxidation of 
four molecules of H2 by several hydrogenases. Four electrons are provided through the reduced form 
of the coenzyme F420, which is regenerated by Frh. Frh is a heterotrimeric enzyme composed of the 
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eLife digest Many microbes grow by producing methane gas from carbon dioxide and 
hydrogen gas, and enzymes known as hydrogenases play important roles in this metabolic process. 
The production of methane in these microbes depends on a nickel-iron hydrogenase called Frh 
adding electrons to a coenzyme called Fi2o- This hydrogenase cleaves a hydrogen molecule Into two 
electrons, which are transferred to the F420 coenzyme, and two protons. The reduced form of F420 is 
then used for several reactions in the methane production process. This process, which is known as 
methanogenesis, provides the microbes with energy. 

Nickel-iron hydrogenases can be divided into five different groups, but researchers have been 
able to determine the detailed structures of the enzymes in just one of these groups. All nickel-Iron 
hydrogenases contain at least two subunlts: a large subunit with a catalytic center composed of 
both nickel and Iron ions and a small subunit that contains three iron-sulfur clusters. Frh — which is 
short for F42o-reduclng nickel-iron hydrogenase — Is known to have a third subunit comprising an 
extra Iron-sulfur cluster and a coenzyme called FAD that allows it to Interact with the F420 coenzyme. 
However, until now, little was known about the detailed structure of the Frh enzyme. 

Mills et al. have used electron cryo-mlcroscopy (cryo-EM) to determine the structure of Frh when 
it Is on Its own, and also when it Is bound to F420. This technique involves freezing a solution of the 
enzyme in a thin layer of ice and recording an Image of this layer in an electron microscope. By 
combining a large number of images, each of which contains many identical enzymes in different 
orientations. It is possible to determine the 3-dlmenslonal structure of the enzyme. 

Mills et al. found that Frh forms a very large tetrahedral complex that contains six Frh dimers. 
And by comparing the structure with and without F420, they Identify a pocket near the FAD coenzyme 
that the F420 coenzyme binds to. They also identify a fold in the third subunit that allows proteins to 
bind both FAD and F420. The work demonstrates the potential of cryo-EM to elucidate structures 
that cannot be determined by other approaches. 
DOI: 10.7554/eLife.0021 8.002 



large subunit FrhA (45 kDa) with a binuclear [NIFe]-center, the small subunit FrhG (26 kDa) with three 
[4Fe4S] clusters and the iron-sulfur flavoprotein FrhB (31 kDa) with a [4Fe4S] cluster and FAD, and the 
F42o-binding site. Frh accounts for about 1% of the cytoplasmic protein and has been shown by electron 
microscopy to form a huge complex of similar appearance in all species investigated: Methanococcus 
voltae {Muth et al., 1987), Methanospirillum hungatei {Sprott et al., 1987), Methanobacterium 
thermoautotrophicum AH {Wackett et al., 1987), and Methanobacterium thermoautotrophicum 
Marburg {Braks et al., 1994). From negative stain Images, the complex was interpreted as a flat cylinder 
consisting of an octamer {Wackett et al., 1987) or hexamer {Braks et al., 1994) of the heterotrimer. 

Recently, developments in instrumentation and image processing software have made it possible 
to determine structures of large macromolecular complexes to near-atomic resolution by cryo-electron 
microscopy. Most successful have been studies of large Icosahedral viruses {Yu et al., 2008, 2011; Liu 
et al., 2010; Wolf et al., 2010; Zhang et al., 2010a; Chen et al., 2011; Settembre et al., 2011), but 
also a few smaller complexes in the 1 MDa size range of lower symmetry have yielded high-resolution 
structures, notably GroEL {Ludtke et al., 2008) and archaeal chaperonins {Cong et al., 2010, 2012; 
Zhang et al., 2010b). We present here the structure of the Frh complex from the hydrogenotrophic 
methanogenic archaeon, Methanothermobacter marburgensis, by cryo-electron microscopy. The Frh 
complex is a 1 .2-MDa dodecamer with tetrahedral symmetry. The location of all cofactors were Identified, 
and the backbone of the three proteins was traced, including FrhB that has a novel fold. 

Results 

The Frh complex was highly purified from M. marburgensis under strict anaerobic conditions in the pres- 
ence of FAD. Its homogeneity was confirmed by electron microscopy with negative staining. We found 
that all projections of the Frh complex are ring-shaped with a diameter of -16 nm {Figure 1A). Class 
averages showed twofold and threefold symmetry {Figure IB), which is consistent with tetrahedral sym- 
metry. We conclude that the complexes consist of six dimers of FrhABG arranged into a shell around a 
hollow core, not In a cylindrical arrangement as in previous EM models {Wackett et al., 1987; 
Braks et al., 1994). The dodecamer has a total MW of 1.2 MDa. Iterative structural refinement of a 
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Figure 1. Cryo-electron microscopy and image processing. (A) A representative area of an electron micrograph taken 
at 200 kV on an FEI Polara. The defocus was determined as 2.05 pm. The scale bar represents 25 nm. (B) Representative 
class averages and corresponding reprojections of the model for the final refinement iteration. The first and last 
images represent a view down the twofold and threefold axis of the tetrahedron, respectively. The scale bar represents 
10 nm. (C) Euler angle distribution for the final reconstruction. Each cylinder represents one class average in the 
asymmetric triangle (1 /1 2 of the tetrahedron); the height of the cylinder is proportional to the number of particles in 
the class. The equal distribution shows that there are no preferred orientations for the Frh complex. 
DOI: 10.7554/eLife.0021 8.003 



84,000-particle cryo-EM dataset, applying tetrahedral symmetry, resulted in a 3D reconstruction {Figure 2) 
at a resolution where the pitch of a-helices and many side chains are recognizable in the map, and 
p-strands are clearly separated [Figure 2E-G, SB). The globular protein complex showed no preferred 
orientations in the ice, resulting in a homogeneous coverage of angles [Figure 1Q. The map shows six 
distinct dimers, all clearly separated from one another [Figure 2C), sitting on the faces of a cube. 
Each monomer contains a series of five high-density features separated by ~1 0 A, forming an arc in the 
interior of the protein [Figure 2D and Video 1). These features are easily assigned to the four Fe-S clus- 
ters and the [NiFe] center; the latter is recognized by its smaller size. Although at first glance the position 
of the FrhABG heterotrimers in the dimer is not obvious, the path of the electron transfer chains provides 
definitive clues to the location of an FrhABG heterotrimer. FrhABG contains one approximately 35-A-long 
helix on the outer surface and a bundle of four up to 40-A-long helices at the interface of different 
dimers, running from the outside to the interior of the complex. Other helices are considerably shorter. 

The large subunit FrhA has considerable homology to the known large subunit structures [Figure 3), 
and is characterized by four long a-helices. Docking the large subunit of the group 1 enzyme into the 
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Figure 2. High-resolution cryo-EM map of Frh. (A) view down the twofold axis. Each of the 12 FrhABG heterotrimers is 
shown in a different color (B) The same view as (A) with the two trimers at the front removed, (C) view down the 
threefold axis. (D) Close-up of the front dimer in (A) shown as a transparent surface with the high densities of the 
metal centers in red. (E) A 10-A thick slice of the map with a ribbon model for one FrhABG subunit (the green one 
in A and D) superimposed. In this and other figures, FrhA is green, FrhG purple, and FrhB slate-blue. At this level, a 
complete chain of cofactors can be seen: the [NiFe] cluster in FrhA (green and brown spheres), three FeS clusters in 
FrhG and one in FrhB (brown and yellow spheres), and the FAD in FrhB (yellow sticks). Details of the map with a 
full-atom model superimposed: (F) Part of the a-helix bundle in FrhA, (G) Helix hairpin in FrhA coordinating the 
[NiFe] center (green and brown balls), (H) Part of FrhG and FrhB showing an FeS cluster and FAD. In (E-H), the map 
was filtered at 4.5-A resolution and sharpened using a B-factor of -54 A^. 
DDI: 10.7554/eLife.0021 8.004 



map provides a good fit for the a-helix bundle, with the [NiFe] center correctly situated in its identified 
location and many other structural elements of FrhA matching closely, including two three-stranded 
antiparallel (S-sheets. FrhA is about 15 kDa smaller than the large subunits of known structure; all 
missing elements are located on the periphery. A structural model for FrhA {Figure 4A) was created 
based on this fit and a sequence alignment with the protein from Desulfovibrio vulgaris Hildenborough 
{Figure 3). All but the first three amino acids in the sequence could be traced in our maps. Densities 
for large side chains confirm the structural assignment {Figures 2, 5 and Video 2). The [NiFe] center 
is located near four cysteine residues, Cys 63, 66, 380, and 383, which are fully conserved in 
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all [NIFe]-hydrogenases. The group 1 [NiFe]- 
hydrogenases contain another ion, typically either 
Mg2+ or Fe^+'^^ {Higuchi et a/., 1997; Garcin et a/., 
1999). In our map, there is strong density at this 
position {Figure 5D), which is part of the structurally 
conserved region of FrhA, and its surrounding 
residues, Glu44, Glu229, and His386, are fully 
conserved in FrhA {Figure 3A), so we interpret 
this density as an ion. One of the helices of the 
four-helix bundle has a pronounced kink. This 
helix was interpreted as residues 166-193, and a 
proline residue (Pro180) not present in the group 
I proteins was found exactly at the position of the 
kink, supporting the assignment. No density is 
present beyond His386, and there is no room in 
the structure for the final 19 amino acids of the 
protein. The C-terminal residues after the equiva- 
lent histidine in [NiFe]-hydrogenases are enzymati- 
cally removed during maturation of the enzymes 
{Menon et a/., 1993; Theodoratou et al., 2005). 
We can conclude that FrhA is posttranslationally 
modified in the same way, likely by FrhD, a 
homologous endopeptidase contained in the Frh 
operon FrhADGB {Liesegang et al., 2010; 
Thauer et al., 2010). 

The fit of the small subunit of the group 1 
enzyme reveals structural homology to the 
N-terminal domain of FrhG including the proximal 
Fe-S cluster. Some peripheral helices present in 
the group 1 protein are missing in the FrhG structure, and the similarity is generally lower than for 
the large subunit {Figure 6C). The N-terminal domain is characterized by a conserved parallel 
P-sheet flanked by short a-helices {Figure 4B). In the group 1 [NiFe]-hydrogenases, the proximal 
Fe-S cluster is located at the C-terminal end of the P-sheet and is coordinated by four cysteines. 
In the primary structure of FrhG, one of the conserved cysteine residues for the coordination of 
this cluster is replaced by Asp60 {Figure 6A). The bacterial hydrogenases as well as most of the 
archaeal species have a cysteine at this position, but aspartates have been identified before as 
ligands of [4Fe4S] clusters {Calzolai et al., 1995; Muraki et al., 2010; Gruner et al., 2011). 
Accordingly, we positioned this aspartate near the proximal Fe-S cluster. 

The C-terminal domain of FrhG has no homology to the group 1 [NiFe]-hydrogenases. Instead, the 
sequence indicates a ferredoxin domain, with two CxxCxxCxxxC sequences coordinating two [4Fe4S] 
clusters {Figure 6A) {Alex et al., 1990). Accordingly, we modeled the C-terminal domain of FrhG on 
ferredoxin {Figures 4B, 5A and Video 3). In our model, the residues 235-244 form a P-hairpin with Gly240 
at the turn. The two P-strands were predicted and the hairpin is clearly visible in the density, confirming the 
assignment. Unlike FrhA, which could be traced over its full length, density of FrhG for a 12-residue stretch 
(188-199) is missing {Figure 6A). This region is located at the surface and may be disordered. FrhG as 
deduced from the M. marburgensis genome sequence is 275 residues long {Liesegang et al., 2010), but 
no density is present for the 45 N-terminal amino acids. The frhG gene most probably starts with the 
initiation codon GTG annotated as Val39 in the genome sequence (Fox et al., 1987; Alex et al., 1990), 
which means that density is missing in the map for only the first seven amino acids of FrhG. 

FrhB contains one Fe-S cluster, an FAD, and the F42o-binding site. There is no homolog of known 
structure. FrhB was therefore traced ab initio, aided by secondary structure predictions {Figure 7) 
and recognizable features in the map that included 10 a-helices, 1 1 P-strands, and the Fe-S cluster. 
The most prominent feature of FrhB is a ~35 A surface helix that corresponds to the predicted 
24-residue C-terminal helix. 

We examined the F42o-binding proteins of known structure in the protein databank, and found no 
structural homology with FrhB. All these proteins have cofactors different from FAD, for example. 




Video 1. The Frh map at different density levels. The 
Frh map contains high densities corresponding to the 
metal clusters. The map, filtered to 5.5 A for clarity, is 
shown at increasing isosurface levels. At density level 1 
(gray), the whole complex is seen; levels 2 (yellow) and 
3 (gold) show decreasing protein density. At level 4 
(orange), only a string of five densities remains, 
corresponding to the four I4Fe4S] clusters and the 
[NiFe] cluster The FeS clusters, but not the INiFe], are 
still visible at level 5 and 6 (red). 
DOI: 10.7554/eLife.0021 8.005 
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Figure 3. Sequence and secondary structure of FrhA and a comparison with homologous proteins. (A) Alignment 
of the [NiFe]-hydrogenase large subunit from Desulfovibrio vulgaris Hildenborough (first line) with 
Methanothermobacter marburgensis FrhA (third line). The second line shows the identical (|) and similar (:) amino 
acids. The alignment was done with EMBOSS needle (http://www.ebi.ac.uk/Tools/psa/emboss_needle/) and 
manually adjusted after fitting the FrhA structure. The fourth line shows the consensus sequence of 19 archaeal 
FrhA species. Identical amino acids in capitals, similar ones in lower case (h: hydrophobic; s: small (GAS); I: large 
(LIFYHW); a: aromatic (FYWH); z: T or S; n: negative, D or E; p: positive, R or K). The a-helices are highlighted in 
green, |3-strands in blue. The [NiFe] ligands are in orange in the ligands of the third ion in red. (B) Comparison 
of the FrhA model (rainbow coloring from blue to red) and the group 1 [NiFe]-hydrogenase large subunit from 
Desulfovibrio vulgaris Hildenborough (gray) (pdb 2wpn) {Marques et al., 2010) in two different orientations (left 
and right). The [NiFe] center and another ion (spheres) overlap as do the structural elements around them, as well 
as two 3-stranded |3-sheets and the lower part of the four-helix bundle. Differences are confined to the periphery. 
DOI: 10.7554/eLife.0021 8.006 
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Figure 4. Models of Frh subunits. (A) FrhA, (B) FrhG, and (C) FrhB in rainbow colors from blue (N-terminus) to red 
(C-terminus). (D) Stereo pair of the FrhABG heterotrimer with the string of four bound Fe-S clusters (orange/ 
yellow), the [NiFe] binuclear center (green/orange) at the top, and the FAD (red) and F420 cofactors (green) below. 
FrhA is green, FrhG is purple, and FrhB is slate-blue. 
DOI: 10.7554/eLife.0021 8.007 



methylene-HiMPT or NADP (Ceh et a/., 2009). Also no homologous FAD-binding proteins were 
found. We conclude that FrhB has a novel fold {Figure 4C and Video 4). The Fe-S cluster is most probably 
a 4Fe-4S cluster because it is surrounded by four cysteines (Cys 104, 134, 192, and 195), which are all 
fully conserved {Figure 7). The core of FrhB is a mixed six-stranded p-sheet. A density located ~8 A from 
the Fe-S cluster could be interpreted as the isoalloxazine ring of the FAD cofactor. The pyrophosphate 
moiety is located at the N-terminal end of an a-helix, which is a common FAD-binding motif where the 
dipole of the helix compensates the negative charge {Dym and Eisenberg, 2001) {Figure 8), and 
nearby density for the adenine moiety is visible. The FAD has a bent conformation, with the adenine 
and isoalloxazine ring only 4-5 A apart {Figure 8) near the highly conserved loop Lys75-Tyr76. The FAD 
is completely surrounded by highly conserved protein regions {Figure 7 and Video 5), further verifying 
the backbone trace. 

To identify the substrate-binding site, we imaged the FrhABG complex in the presence of a 
large excess of F420, and calculated maps, which were refined to a similar resolution. At all refinement 
stages, the map showed an additional density in a region that was empty in maps of Frh without 
substrate {Figure 8A,B). This density runs approximately parallel to the FAD isoalloxazine ring, at 
a distance of ~4 A. We conclude that it represents the isoalloxazine ring of the F420 substrate. It is 
seated in a large pocket surrounded by highly conserved residues, with easy access to the surface 
of the complex. The F420 density is close to the small residues S209 and V210 in a loop between 
two p-strands, which are part of the longest conserved region of FrhB {Figure 7). A predicted 
three-stranded P-sheet (148-172) flanked by two short helices {Figure 7) fits nicely in density. The 
only conserved region of this sheet, 163IGKGK, forms a turn close to the position of the isoalloxazine 
ring of F42o. The location between this ring and the surface of the complex suggests that one or 
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Figure 5. Details of the model and map. (A) The ferredoxin domain of FrhG. The eight cysteine residues surrounding 
the two [Fe4S4] clusters are shown in green. (B) P-Sheet in FrhB. The strands are clearly separated. (C) The dimer 
interface of FrhA inside the particle with two salt bridges Asp312-Arg314. The two FrhA molecules are viewed from 
inside the complex. (D) The C-terminal three-stranded P-sheet in FrhA with the [NiFe] center (green and brown) and 
another ion (gray) near the C-terminal His386. The conserved residues Glu44, Glu229, and His386 coordinate the ion. 
DOI: 10.7554/eLife.0021 8.008 



both of the conserved lysines may interact with the phosphate group of F42o- The long a-helix on 
the surface that was interpreted as the C-terminal helix of FrhB is not very well resolved. The residues 
in the helix are mostly not conserved in FrhB, except at the N-terminal end {Figure 7). In our 
model, the conserved residues are located near the loop 164GKGK mentioned above. This region 
also contains several other conserved residues {Video 5} and likely constitutes the access channel 
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for the substrate F42o- The map of Frh in the 
presence of substrate is similar to the substrate- 
free map, indicating that F420 binding does not 
involve major conformational changes. 

The residues coordinating the Fe-S cluster, FAD, 
and F420 are conserved in the FrhB family, which 
includes subunits of F42o-dependent sulfite reduct- 
ase (Fsr), F42oH2:quinone oxidoreductase (FqoF), 
F42oH2:phenazine oxidoreductase (FpoF), F420- 
dependent glutamate synthase and formate dehy- 
drogenase (FdhB) {Johnson and Mukhopadhyay, 
200Sj {Figure 9). Some of these proteins contain 
additional domains, like an N-terminal ferredoxin 
domain (which in Frh is part of FrhG), or a C-terminal 
domain containing a binding site for the substrate 



Video 2. FrhA. The model of FrhA is shown of the electron transfer from F420. Clearly, these 

superimposed on the map. enzymes share a common fold with FrhB. 
DOI: 10.7554/eLife.0021 8.009 

Discussion 

Recent advances in instrumentation and image 
processing procedures have made cryo-electron microscopy of isolated macromolecular complexes a 
powerful technique in structural biology. After the first virus structure at subnanometer resolution 
{Bottcher et a/., 1997) showed the feasibility of the method, icosahedral viruses were reconstructed 
in recent years to better than 4 A resolution (Yu et a/., 2008, 2011, Zhang et a/., 2008, 2010a; 
Liu et al., 2010; Wolf et al., 2010; Chen et al., 2011; Settembre eta/., 2011; ), taking advantage 
of the 60-fold symmetry and huge size (~20 to 1 50 MDa), which together contribute to the relatively 
high signal-to-noise ratio for these particles that makes the determination of the orientation 
unambiguous. For smaller complexes, the techniques have been pioneered with the 800-kDa 
GroEL complex with D7 symmetry, reaching subnanometer resolution for the first time in 2004 
(Ludtke et al., 2004) and ~4 A in 2008 {Ludtke et al., 2008). Subsequent reconstructions to this 
resolution were achieved for archaeal {Zhang et al., 2010b) and mammalian {Cong et al., 2010) 
chaperonins. The chaperonins have different conformational states, and usually only one of these 
gave a high-resolution map, whereas for the others the resolution was limited due to conformational 
flexibility. The appearance of our reconstructions is similar to the chaperonins; in agreement with 
this, Fourier shell correlation indicates a resolution of 3.9 and 4.0 A for the two maps, respectively 
{Figure 10A). Recently it has been noted that the resolution of such reconstructions are usually 
overestimated due to correlation of noise at high resolution {Scheres and Chen, 2012). Using a 
'gold standard' procedure to estimate resolution, where two half-data sets are independently refined 
to a reference not containing high-resolution information (see 'Materials and methods'), yielded 5.5 
A {Figure 10B). However, an objective indication of (local) resolution is provided by the appearance 
of secondary structure: the a-helical pitch of 5.4 A is very well defined in many regions {Figure 2F,G), 
as is the separation of P-strands (axial distance -4.8 A) {Figure SB). Comparison of these areas 
with maps calculated from the model low-pass filtered to different resolutions suggest a resolution 
-4.5 A. The overall resolution estimate is an average of these well-resolved protein regions and 
other regions that are less well resolved, especially those at the periphery of the complex. An FSC 
curve between the experimental map and a map calculated from the model indicated a resolution of 
5.8 A {Figure 10Q; this figure also constitutes an average over the whole structure, with the model 
probably best for internal a-helices and worst for loops at the periphery. Since the model was 
derived from the map, and not from an independent experiment, this measure does not properly 
reflect the resolution of the map, but rather describes how well our admittedly incomplete model 
explains the map densities. It should be noted that we used tight protein structure restraints during 
modeling, in order to prevent overfitting to noise while creating an unrealistic model that would 
have yielded a higher model-to-map correlation. 

Several factors contributed to the successful structure determination of Frh. The tetrahedral symmetry 
means that there are 12 asymmetric units in each particle, less than the 14 and 16 in the D7 and D8 
symmetric chaperonins, but the resulting ball-like shape ensures that the complex has no preferential 
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Figure 6. Sequence and secondary structure of FrhG and a comparison with homologous proteins. (A) Alignment 
of the [NiFej-hydrogenase small subunit from Desulfovibrio vulgaris Hildenborough (first line) with 
Methanothermobacter marburgens/s FrhG (third line) and ferredoxin from Peptostreptococcus asaccharolyticus 
(pdb 1 dur) (last line). The second and fourth lines show the identical (|) and similar (:) amino acids. The alignment 
was done with ClustalW. The a-helices are highlighted in green, P-strands in blue. Ligands of the proximal [4Fe4S] 
cluster are in yellow, of the medial cluster in orange, and of the distal cluster in red. Amino acids not seen in the Frh 
structure are in gray font. The first 38 amino acids of the FrhG sequence deduced from the M. marburgensis 
genome are probably not part of the protein (see main text) and are shown in light gray font. (B) Consensus 
sequence of archaeal FrhG species. Capitals: conserved residues; lower case: similar residues (h: hydrophobic; s: 
small (GAS); I: large (LIFYHW); a: aromatic (FYWH); z: T or S; n: D or E; p: R or K). (C) Comparison of the FrhG model 
(rainbow coloring from blue to red) and the [NiFe]-hydrogenase small subunit from Desulfovibrio vulgaris 
Hildenborough (pdb 2wpn) (gray) in two different orientations (left and right). The structure near the proximal 
[4Fe4S] cluster (blue-green) is conserved, but the periphery diverges. The ferredoxin domain (yellow to red) 
containing the medial and the distal iron-sulfur clusters is not homologous and the clusters do not overlap. 
DOI: 10.7554/eLife.0021 8.010 
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orientation in the ice {Figure 1D), unlike the 
cylinder-shaped chaperonins that tend to be 
oriented with their symmetry axis parallel or 
perpendicular to the ice {Ludtke et al., 2001; 
Zhang et al., 2011). From the quality of the 
reconstructions, we conclude that the complex 
is very rigid. The maps of Frh with and without 
substrate are of similar resolution and quality, 
implying that any protein conformational changes 
due to F420 binding are too small to detect at 
the present resolution. The density for F420 is 
weak, but its position was inferred at an early 
stage of the analysis from the fact that maps of 
Frh without the substrate were consistently 
empty and maps with F420 always showed some 
density here. The later tracing of the amino 
acid chain of FrhB and the FAD confirmed the 
assignment: F420 is in van der Waals distance 
from its hydride donor, the FAD isoalloxazine 
ring, and the residues in the vicinity are highly 
conserved In FrhB and in the FrhB family of 
F42o-binding proteins {Figures 7, 9 and Video 5). 
While the chain trace for FrhA and FrhG was relatively straightforward using the group I hydro- 
genase homolog, FrhB was built ab initio. The visibility of the FeS cluster density {Figure 2D and 
Video 1) made it possible to use the four Fe-S ligand cysteines together with a secondary struc- 
ture prediction as an initial guide; large side chain densities helped to establish the correct 
sequence assignment. The secondary structure prediction was fully confirmed, and weak or dis- 
continuous densities in the map were almost invariably found to represent glycine residues, further 
indicating a correct assignment. Density for the FAD cofactor only became clear when most of the 
protein chain was placed; density for the ribose and ribltol are missing {Figure 8). The pyrophos- 
phate group is closely associated with the N-terminal a-helix 26-64 of FrhB. The N-terminus of this 
helix is one of the most conserved regions in the FrhB family {Figures 7, 9), and its role in com- 
pensating the negative charge of the phosphates is a common motif in FAD-blnding proteins 
(Dym and Eisenberg, 2001). Conserved regions of FrhB were found near the FAD and surrounding 
the F42o-binding site and its putative entrance channel {Video 5). Of the three proteins, only FrhG 



Video 3. FrhG. The model of FrhG is shown 
superimposed on the map. 
DO!: 10.7554/eLife.0021 8.011 
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Figure 7. Sequence and secondary structure of FrhB. Light gray font: No density; green highlight: a-helix; gray 
highlight: P-strand. Second line: Consensus sequence of FrhB species. Capitals: conserved residues; lower case: 
similar residues (h: hydrophobic; s: small [GAS]; I: large [LIFYHW]; a: aromatic [FYWH]; z: T or S; n: negative [D or E]; 
p: positive [R or K]). In the consensus sequence, the residues for coordination of the iron-sulfur cluster, FAD, and 
F420 are highlighted in orange, cyan, and yellow, respectively. 
DOI: 10.7554/eLife.00218.012 
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is missing density for an internal stretch, resi- 
dues 188-199 {Figure 6A). The flanking resi- 
dues of this stretch are located less than 10 A 
apart on the surface of the complex and they 
probably form a flexible loop. FrhA and FrhB 
were traced over their full length. All three 
proteins follow the rule that hydrophobic resi- 
dues form the core of the protein and charged 
residues face the outside. We found several 
pockets lined with hydrophobic side chains, not 
just in FrhA and FrhG, where these features 
are shared with the bacterial protein used as 
a template for modeling, but also in FrhB. 
FrhB also contains a highly hydrophobic helix 
(27GIVTGLLAYAL), which is buried completely 
inside the subunit. These features all indicate a 
correct chain trace. 

Our structure shows that each heterotrimer 
encompasses a complete electron transfer chain 
that runs from the [NiFe] cluster in FrhA via the three Fe-S clusters in FrhG and one Fe-S cluster 
and FAD in FrhB to F420 {Figures 4D, 11 and Video 6). The cofactors describe a curvilinear path in 
the protein interior, almost parallel to the surface of the complex, with distances of ~1 0 A between 
neighbors, except FAD and F420, which are in van der Waals contact to enable hydride transfer 
{Figure 8C). The predicted F42o-binding site is located on the outside of the complex, which facili- 
tates access of the cytosolic F420. The two electron transfer chains in a dimer are separated by at 
least 17 A and are thus independent of each other, indicating that the functional unit is one 




Figure 8. Electron-transfer chain and F42o- Maps in the absence (A) and presence (B) of the substrate F420 differ in a 
region near a conserved loop between two P-strands near the FAD (carbons in yellow). The isoalloxazine ring of F420 
fitted into this density is shown in pink. (C) Electron transfer chain with minimal distances indicated. 
DOI: 10.7554/eLife.00218.014 




Video 4. FrhB. The model of FrhB is shown 
superimposed on the map. 
DOI: 10.7554/eLife.00218.013 
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FrhABG heterotrimer and there Is no coopera- 
tivity in the complex. The three proteins in the 
Frh heterotrimer form a tight complex with 
extensive subunit interfaces and a large number 
of subunit interactions, including two salt 
bridges between loops of the two FrhA mole- 
cules in a dimer, from Asp312 of one FrhA and 
Arg314 of the other {Figure 5C). The hetero- 
trimer is quite slender compared to the two- 
subunit [NiFe] hydrogenases of known structure, 
with a small contact surface especially between 
FrhB and the other two subunits {Figure 12), 
and complex formation may be necessary for 
stability of the complex and keeping the metal 
clusters at optimal distances. The same quater- 
nary structure of Frh has been observed not 
only in the closely related thermophile M. 
thermautotrophicus {Wacl<ett et al., 1987) but 
also in the phylogenetically distant mesophilic 
M. vo/tae {Muth et al., 1987) and M. hungatei 
{Sprott et al., 1987), so it is clear that the for- 
mation of a large complex plays an essential 
role for the function of Frh. 

Materials and Methods 
Cultivation 

Methanothermobacter marburgensis (DSM 2133) was obtained from the Deutsche Sammlung von 
Mikroorganismen (DSMZ, Braunschweig, Germany). The archaeon was grown anaerobically at 65°C on 
80% H2/20% COa/O.IX H2S in a 12-L fermenter containing 10 L complete mineral salt medium 
{Schonheit et al., 1980). Cells were harvested by the use of a continuous-flow centrifuge under 
anoxic conditions at the late exponential phase and stored at -80°C. 

Purification of the Frh complex from M. marburgensis 

Purification was performed under strictly anaerobic conditions at 18°C in an anaerobic camper 
(Coy Laboratory Products, Grass Lake, Ml). All buffers used contained 2 mM DTT and 25 |jM FAD. 
Cell extracts were routinely prepared from 20 g (wet mass) of M. marburgensis cells. The cells 
were suspended in 35 ml 50 mM Tris/HCI pH 7.6 (buffer A), and passed four times through a 
French pressure cell at 125 MPa. Cell debris were removed by centrifugation at 15,000xg for 
30 min. The supernatant, designated cell extract and containing -1000 mg protein, was adjusted 
to 800 mM ammonium sulfate in buffer A and stirred for 20 min. The cell extract was applied to a 
Phenyl Sepharose 6 Fast Flow column (6x10 cm) equilibrated with 800 mM (NH4)2S04 in buffer A. 
Protein was eluted by a (NH4)2S04 step gradient in buffer A: 800 mM (NH4)2S04 for 250 ml, 
200 mM (NH4)2S04 for 250 ml, and 0 mM (NH4)2S04 for 50 ml (flow rate: 8 ml/min). The Frh activity 
was eluted in the 0 mM (NH4)2S04 fractions. The protein solution was concentrated with ultrafiltration 
by Amicon filters (100-kDa cutoff) to 4-5 ml, which were then applied to a Sephacryl S-400 HR 
column (2.6 x 60 cm) equilibrated with buffer B (buffer A -1- 150 mM NaCI). The Frh activity was 
eluted after washing the column with 180 ml buffer B (flow rate: 1 ml/min). The pooled fractions 
were concentrated with ultrafiltration by Amicon filters (100-kDa cutoff) to 20 mg/ml. CHAPS was 
added to the concentrate (48 mM final concentration), and the solution was incubated for 12 hr at 
room temperature with slow stirring. The protein solution was washed 5 times on Amicon filters 
(100-kDa cutoff) with buffer A with 4 mM CHAPS (buffer A2) and applied to a MonoQ column 
(1x8 cm) equilibrated with buffer A2. Protein was eluted by a NaCI linear gradient in buffer A2: 
0-400 mM NaCI in 25 ml and then 400-600 mM NaCI in 20 ml (flow rate: 0.8 ml/min). In the linear 
gradient of 400-600 mM NaCI, the active Frh was eluted from the column at a concentration of 
540 mM NaCI. The protein solution was concentrated on Amicon filters (100-kDa cutoff) to 2 ml 




Video 5. Conserved residues in FrhB. The FrhB model 
is shown with conserved residues in red, similar residues 
in orange and unconserved residues in gray (see Figure 7 
for sequence). FAD is yellow and F420 green. Highly 
conserved regions surround the FeS cluster and the FAD 
as well as the F420 isoalloxazine ring. An entrance pathway 
for F420 is suggested as well. 
DOI: 10.7554/eLife.00218.015 
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Cons 
FrhB 
GS 

FqoF 114 
FpoF 74 
Fsr 66 
FdhB 1 



1 

268 



1-clhhAps sQDGGhhzsll-Yhlcc — hnshhhh acs — hhh — s — h 6 8 

MVLGTYKEIVSARSTDREIQKLAQDGGIVTGLLAYALDEGIIEGAWAGPG-EEFWKPQPMVAMSSDEL 62 

-EPLGEYTEILSARA PMFRGQDGGWTALLTYALREGI VDGALVVDRDPAMPWKPVPVLAEDPEDV 3 32 

-NGLGEYIEVVAARSKR FVGQDGAMVTEFTASALEMGI lERAIFVARD — SNWRTRVVTIKTPEQL 177 

-NELANVRKFFAARS KENAGSQDGGVTSGILKSLFKQGKIDCAVGITR — DEKWESKVVLLTSAEDV 137 

-IREKFYEKYYYAKS DIEGQDGGWTAFLKYLLENGKIDGAIVVG DECWKPVSLVVQNAEDL 126 

MKYLLARATDEEIQRKGECGGAVTAIFKYMLDKEVVDAVLTLERG-YDVYDGIPVLLEDSSGI 6 2 



Cons — h — zKY h — L phshhh-PCQ — shp-hc 

FrhB 6 9 KAAAGTKYTFSPNVMMLKKAVRQYGIEKLGTVAIPCQTMGIRKMQTYPFGVRF LADKIK 127 

GS 333 VRAAGTKYSVCPILKVLKE351 IMSRYAMVGTPCQITAATLMKEYN GEFPVE 29 

FqoF 178 YDRKITGTKYSYADVLPALKEAVLKSEAVGFVGTPCMVSAVRKMQQAFKKFER VK 2 32 

FpoF 138 EKVRGTKYTSDPVVAALREAFEKY — DRIAVVGVPCQAHSARLIREN VSEKIV 188 

Fsr 127 LKTAKSKYAIST-LDALRKAGEMG-LEKVAVVGLPCQINGLRKLQYFPYHAKHDLELGRNGKPVKLPKIE 194 

FdhB 63 ESTCGSLHCAPTMFGDLISRYLSD — MRLAVAVKPCDAMAIRELEKRHQ IDPDKV 115 

Cons h-hGllC 1 1 h-hc-V cl-pGpllh-h-n h-lc C 

FrhB 128 LLVGIYCMENFPY-TSLQTFICEKLGVSMELVEKMDIGKGKFWVYTQD-DVLTLPLKETHG YEQAGC 192 

GS 30 LRIGLFCMENFSY-TYLREL-AEAEGVDLRDVSECRIEKGRLWFHLNDGSTVSIPLERARS AMRKNC 94 

FqoF 233 LAIGLFCTENFYH-HDLYKFLLEKANADLRNAVKTDIKKGKFIVEMKDGSKVRIPVKDFEE IIPSGC 295 

FpoF 189 LIIGLLCMESFHHDVMLDKIIPEIMKVKIEDVRKMEFTKGKFWVYTSDGEVHSVPIKDVAK YARNPC 255 

Fsr 195 YLIGLFCTEKFRYDNMKEVLS — KHGIDIEKVEKFDIKKGKLLVYVN-GEKKEFDLKEFEI CSGC 256 

FdhB 116 YKIGLNCGGTLAPVSAREMIETF-YEIDPDDVVSEEIDRGKFIVELRDGSHREISIDYLEEEGFGRRENC 184 

Cons — C — 1-s — sDhshG-hGs-DGaz-h — pTn-Gp-l-c-h hcKh 

FrhB 193 KICKDYVAELADVSTGSVGSPDGWSTVITRTDAGDSIFKQAVEAGLFETKPIEEVKPGLGLLEKLAAQKK 262 
GS 95 SVCMDFTSEQSDVSVGSVGSPQGWSTLIIRTERGRELVDGAAKAGYIETAPITGK — GLKLLEKLASGKK 162 

FqoF 296 KVCQDFAAVESDVSVGSVGSPNRFSTVMVRTEVAKQILDYIREKDYAEFG EAKLDIVQKLCDHKM 363 

FpoF 256 HHCCDYTSVFADISVGSVGAPDGWNSVFIRTDAGEEYFEMVREE MEIMEDPKPGLELVKKLIDMKR 321 

Fsr 257 KMCRDFDAEMADVSVGCVGSPDGYSTIIIRTEKGEEIKNAVELK EGVNLEEIEKLRQLKL 316 

FdhB 185 QRCEIMVPRNADLACGNWGADDGWTFIEVNTERGQEIIEGARSSGYIE AREPSEKMVKIREKIE 248 

Cons 

FrhB 263 EKAEKNIAARKEMGLPTPF 281 

GS 163 EENLEEIQRRESVARPVLYWRVMPGDLYPEEIKDCQFDDLRADVIDVGACVLCGACEASCPEGIVRIND- 231 

FqoF 364 KIHPWPPKKKEEKSEE 379 

FpoF 322 KNNAEHFKEVCKEFSFETGIRDETV 34 6 

Fsr 317 KRFKKEVERRRENNEYVSFYWTADYGGIGKRADGTYFIRVRAKPGGWYKPEEIKEILDIAEEYNAKIKV- 385 

FdhB 249 NAMISMARKFQDKYLDEEYPSLDEWDEYWKRCINCFACRDACPVCFCRECELEKDYLLESDEKAPDPLT- 317 



Figure 9. Alignment of the FrhB family of F42o-binding proteins. FrhB: FrhB of Methanothermobacter marburgensis 
DSM 2133; GS: F42o-dependent glutamate synthase from M. marburgensis, ADL58239; FqoF: F42oH2:quinone 
oxidoreductase subunit F from Archaeog/obus fulgidus DSM 4304, NP_070660; FpoF: F42oH2:phenazine 
oxidoreductase subunit F from Methanosarcina barker! sir. Fusaro, YP_303819; Fsr: N-terminal domain of Fsr, 
F42o-reducing sulfite reductase from Methanocaldococcus jannaschii DSM 2661, Y870_METJA; FdhB: beta 
subunit of formate dehydrogenase from M. marburgensis, YP_003850414. Cons: the consensus sequence of the 
FrhB family. Capitals: conserved residues; lower case: similar residues (h: hydrophobic; s: small [GAS]; I: large 
[LIFYHW]; a: aromatic [FYWH]; z: T or S; n: negative [D or E], p: positive [R or K], c: charged [D, E, R, K]). In the 
consensus sequence, residues indentified in the FrhB structure for coordination of the iron-sulfur cluster, FAD, 
and F420 are highlighted in orange, cyan, and yellow;, respectively. Secondary structure prediction for each of the 
proteins was done with PSIPRED (http://bioinf.cs.ucl.ac.uk/psipred/): green highlight: a-helix; gray highlight: 
P-strand. Highly conserved amino acids are indicated in bold font. Amino acid numbers are shown in red; note 
that the C-terminal 83 residues of glutamate synthase (GS) align with the N-terminus of FrhB and the N-terminal 
182 residues with the C-terminus of FrhB. 
DOI: 10.7554/eLife.00218.016 



and further purified on a Sephacryl S-400 HR column in order to remove remaining smaller particles 
and bigger aggregates. 

F420 was isolated from M. marburgensis by established methods {Shima and Thauer, 2001). 

Cryo-electron microscopy 

3 |jl of a 0.7 mg/ml Frh sample was applied to freshly glow discharged Quantifoil R1/4 grids 
(Quantifoil Micro Tools, Jena, Germany) that had been pretreated for 15 s in chloroform. The 
Igrids were blotted in an FEI Vitrobot using a 2.5 s blotting time at 70% humidity and 10°C and 
plunge-frozen in liquid ethane. For collection of data of Frh in the presence of substrate, 0.5, 1 .0, 
or 10 mM F420 was added to the sample under oxygen-free conditions just prior to grid 
preparation. 
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Figure 10. Resolution estimations and B-factor determination. (A) Fourier shell correlation plot, showing the 
resolution for the substrate-free map (black) and the F42o-containing map (red). At the 0.5 FSC criterion, the 
resolution is 3.9 and 4.0 A, respectively. (B) Gold standard FSC between two half-data sets independently refined 
from a low-resolution model. The resolution at the 0.143 FSC criterion is 5.5 A. (C) FSC between the map and 
model. At 0.5 FSC, the indicated resolution is 5.8 A. (D) Plots of the natural logarithm of the spherically averaged 
structure factor amplitude as a function of the resolution (A"^). Blue: experimental amplitudes; red: Cef-weighted 
amplitudes; green: amplitudes sharpened with a B-factor of -54 A^, as determined by the program embfactor 
[Fernandez et al., 2008) for the resolution zone 10^.5 A. The purple line shows the average scattering amplitude 
for the Frh complex (VNa,„,^J. 
DOI: 10.7554/eLife.00218.017 



Images were collected at liquid nitrogen temperature on an FEI Tecnai Polara operated at 200 kV. 
Before images were recorded, the microscope was carefully aligned in an iterative process to correct 
for objective astigmatism and beam tilt by coma-free alignment (G/aeser et a/., 2011). The corrections 
were carried out at a dose of 15 e /A^ and at half the defocus value used for collecting the images, 
and were repeated for each grid square from which images were collected. The alignments were done 
with a Gatan 4k x 4k CCD using unbinned images. 

Images were recorded on Kodak SO-163 film at a magnification of 59,000x with a dose of 
10-15 e /A^ at a defocus of 1 .5-2.8 |jm. The film was developed for 12 min in full-strength Kodak D-19 
developer and fixed for 8 min in Kodak Rapid Fix. Five hundred and six negatives were collected of 
substrate-free Frh and 512 of Frh with Fn2o- Films that showed obvious flaws (too thin ice with particles 
just around the edge of the Quantifoil hole, too high particle concentration, broken ice, obvious 
drift) were discarded. Four hundred and two negatives of substrate-free Frh were scanned and 277 of 
F420 with Frh. The latter preparation had a higher protein concentration, resulting in many films with 
too many overlapping particles. This was more than compensated by the higher number of parti- 
cles per film. 

Image processing 

Films were digitized on a Zeiss Photoscan scanner with a pixel size of 7 jjm, corresponding to 1.14 
A on the specimen as calibrated with fatty acid synthase (Gi'pson et al., 2010). Particle selection 
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Figure 11. Model of the dodecameric Frh complex. Each FrhABG heterotrimer is colored differently with the same color scheme as the map in Figure 1. 
(A) View down the twofold axis, (B) view of the FrhB trimer, and (C) view of the closely packed FrhA trimer as Figure IC. 
DOI: 10.7554/eLife.00218.018 



was done semlautomatically with the Boxer module from EMAN {Ludtke et a/., 1999) and data 
processing with EMAN2 (Tang et a/., 2007). The contrast transfer function (CTF) of the selected 
particles from each film was determined with EMAN2 and images that showed visible Thon rings 
in the power spectrum to high resolution and no indication of drift or astigmatism were selected 
for further processing and their CTF was corrected by phase flipping. A tetrahedral starting model 
was created by the EMAN2 program e2initialmoclel from a number of class averages. This model 
was iteratively refined using the main EMAN2 program eZrefine, which determines the 3D orienta- 
tion of each particle by comparison to a set of projections of the current 3D reference map. 
Particles in the same orientation are aligned and averaged, and a new 3D map is constructed from 
the averages that are then reprojected to create references for the next refinement cycle {Ludtke 
et a/., 1999; Tang et a/., 2007). Tetrahedral (T) symmetry was applied throughout. In the initial 
refinement steps, the data were binned to a pixel size of 2.28 A. After the resolution reached -10 A, 
the unbinned data were used. At every refine- 
ment step, the 30% worst members of each class 



were discarded. Because of the high symmetry 
and globular shape of the particle, all class aver- 
ages had similar quality (see Figure IB) and were 
used for each reconstruction. The final data set of 
substrate-free Frh contained 84,000 particles 
from 101 negatives of the 402 scanned nega- 
tives. Between iterations, the projection angle 
was varied. In the last refinement step an angle of 
0.9° between consecutive reference reprojections 
was used, yielding 2134 reference images. 

The F420 data were refined with a high-resolution 
map of substrate-free Frh as a starting refer- 
ence map. Ninety seven thousand particles were 
selected from 80 negatives. Smaller datasets of 
Frh with 0.5, 1.0, and 10 mM F420 were used ini- 
tially, but no significant differences were seen 




and the datasets were merged. 

The resolution of the maps was estimated Video 6. The Frh dodecamer. A model of the 

using the command eotest in EMAN2, which cal- tetrameric Frh complex, 

culatesthe Fourier shell correlation (FSC) between DOI: 10.7554/eLife.00218.019 
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Figure 12. A comparison of (A) Frh and (B) the group I [NiFe] hydrogenase from Desulfovibrio gigas (2wpn) 
[Marques et a/., 2010) shows a lower intersubunit contact area in the former The large subunit/FrhA is shown in 
green, the small subunit/FrhG in magenta, and FrhB in blue. The left panel shows Frh in an orientation along the 
twofold axis of the complex, as seen from the outside, and the right panel a view 90° rotated, as seen from the 
dimer partner 

DOI: 10.7554/eLife.0021 8.020 



reconstructions made by half-data sets from the odd and even numbered particles. This indicated 3.9 
and 4.0 A for the maps without and with F420, respectively, using the 0.5 FSC criterion {Figure 10A). 
The resolution estimate can be inflated by the overfitting of noise at high resolution {Grigorieff, 2000; 
Scheres and Chen, 2012). To test for this effect, we used the procedure e2refine-evenoclcl in EMAN2, 
where the phase of the starting model are randomized from a cutoff resolution lower than the expected 
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resolution of the map, and then two half-data sets are refined to convergence against this model. 
High-resolution noise will be uncorrelated between the two maps. Subsequently, the FSC between 
the two resulting maps is calculated, and the resolution from these two independent half-data sets 
determined at 0.143 FSC {Rosenthal and Henderson, 2003). This procedure yielded a resolution of 
5.5 A [Figure 10B). B-factors of the reconstructions were estimated with the program embfactor 
{Rosenthal and Henderson, 2003; Fernandez et al., 2008) in the resolution range 10-4.5 A as 54 A^ 
for the substrate-free map and 60 A^ for the map with F420 (see Figure 10D). 

Model building 

Maps were visualized in Chimera {Pettersen et al., 2004). Homologous protein structures were 
fitted manually in the map followed by a rigid body fit (pdb 2wpn, Marques et al., 2010, for FrhA 
and the C-terminal domain of FrhG; pdb Idur for the ferredoxin domain of FrhG). Models were 
then mutated to the correct amino acid and manually rebuilt in Coot {Emsley and Cowtan, 2004). 
Regions for which no homolog was available, including all of FrhB, was built ab initio in Coot, 
based on secondary structure predictions done using the PSIPRED server {Bryson et al., 200S). 
The conserved cysteines serving as ligands for the Fe-S clusters and the visibility of the Fe-S clus- 
ters in the map gave clear initial hints to the identity of secondary structure elements. Helices were 
built using the command Place helix. Loops and (B-sheets were built by first placing Ca atoms using 
the Baton build module. Full-atom models were built based on the visibility of side chains, which 
led in all cases to an unambiguous assignment. The models were refined using the Regularize 
Zone module. Where possible, side chains were fitted to available density (overall about 55% of 
side chains were visible, including most aromatic residues and arginines). During refinement, tor- 
sion angle, planar peptide, and Ramachandran restraints were used in order not to create a well- 
fitting but unrealistic model. The final model contains 878 residues out of a possible 903, a [NiFe] 
cluster, 4 [4Fe4S] clusters, an FAD and part of an F420 molecule; 88.1% of residues lie in the most 
favored regions of a Ramachandran plot, 8.0% in generously allowed regions, and 3.8% are out- 
liers. A full dodecamer model was generated in Chimera {Pettersen et al., 2004) and a map was 
calculated from this model. An FSC curve between this map and the experimental map indicated 
an overall agreement of 5.8 A at 0.5 FSC {Figure IOC). 

Figures were made using Chimera {Pettersen et al., 2004) and the PyMOL Molecular Graphics 
System {DeLano, 2002). 
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